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Abstract 

In this paper, we study the existence of higher order Poisson type systems. In detail, we 
prove a Residue type phenomenon for the fundamental solution of Laplacian in M", n > 3. 
This is analogous to the Residue theorem for the Cauchy kernel in C. With the aid 
of the Residue type formula for the fundamental solution, we derive the higher order 
derivative formula for the Newtonian potential and obtain its appropriate C'^'" estimates. 
The existence of solutions to higher order Poisson type nonlinear systems is concluded as 
an application of the fixed point theorem. 



1 Introduction and background 

We study the existence of solutions u = {ui, . . . ,un) to the following nonlinear system in M", 
n > 3: 

A"^m(x) =a(a;,u,Vu,...,V^'"M). (1) 

Here m G V-'m represents all j-th order partial derivatives of all the components of u, and 
a := (ai, . . . , ttiv) is a vector- valued function on x and the derivatives of u up to order 2m. Label 
the variables of a by (p_i,po,Pi, • • • ,P2m), with p_i representing the position of the variable x 
and pj representing the position of V%, < j < m. 

The solvability has been one of the central problems in the theory of partial differential 
equations and has been explored widely since the counterexample of Hans Lewy |Lw] in 1957. 
See |Hol] . |NTj . [Mo] . |Ho2] . |BF] . |Ln] . |De] and the references therein. Unlike linear equations, 
there is in general no systematic theory about the solvability for nonlinear equations, not to 
mention nonlinear system. Most recently, [Panlj investigated the existence problem to ([T]) 
through Cauchy Riemann operator in the case when n = 2. In a consequent paper [Pan2j . Pan 
studied the solvability when m = 1 for general n > 3. 

Our main theorems are stated as follows. 
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Theorem 1.1. Let a G C^''^{0 < a < 1). For any given appropriate initial data {cj}o<j<2m-i; 
there exist infinitely many solutions in the class of C^™'" to the initial value system 

A™u(x) = a{x, M, Vm, . . . , V^'^-^m); 
m(0) = Co; 

Vu(0) = ci; (2) 

in some small neighborhood of 0. Moreover, all those solutions are of vanishing order at most 
2m and not radially symmetric. 

Here we call {cj}o<j<2m-i cin appropriate initial data if it satisfies the appropriate 
symmetry conditions as derivatives of any vector- valued function. This is apparently necessary 
for the existence of solutions to the system ([2]). A function u in the class of C'' is said to be of 
vanishing order m (m < k) at if V%(0) = for all < j < m — 1 and V'"u(0) ^ 0. 

We point out, since the solutions obtained in Theorem 11.11 are of vanishing order at most 
2m, they are never trivial solutions. Moreover, since the solutions are not radially symmetric, 
they are not obtained by possibly reducing the system into an ODE system with respect only 
to the radial variable r = 

Due to the flexibility of a. Theorem 11.11 can be used to construct local m-harmonic maps 
from Euclidean space to any given Riemannian manifold. The resulting image in the target 
manifold can be either smooth or singular, depending on the given initial data. 

When a is dependent also on p2m variable, we obtain the following existence theorem with 
some additional assumption on a. 

Theorem 1.2. If a G and a(0) = Vp2„a(0) = Vp2„a(0) = 0, then there exist infinitely many 
solutions in the class of C^™'" (0 < a < 1) to the system 

A'^uix) = a{x, M, Vu, . . . , V^'^u) 

in some small neighborhood of 0. Moreover, all those solutions are of vanishing order 2m and 
not radially symmetric. 

On the other hand, if the system ([1]) is autonomous, i.e., independent of the variable x, then 
there exists solutions over large domains in the following sense. 

Theorem 1.3. If a E and a(0) = Va(0) = 0, then for any R > 0, there exist infinitely 
many solutions in the class of C"^"^'" to the autonomous system 

A'^u = a{u, Vm, . . . , V^'^u) 

in {x G : \x\ < R}. Moreover, all those solutions are of vanishing order 2m and are not 
radially symmetric. 
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We would like to point out, even though the autonomous system in Theorem 11.31 is itself 
translation invariant, none of the solutions obtained there is a trivial translation of the radial 
solution, from the proof of Theorem ll.3[ On the other hand, the regularity of a in Theorem 
11.31 can be reduced to 0^'°" if a is in addition independent of V^^m variable. This will be seen 
from the proof of Theorem 11.21 and 11.31 in section 9 and 10. The fact will be used in some of 
the following examples. 

We also note the neighborhood in Theorem 11.11 where the solutions exist is necessarily small, 
as indicated by the following example of Osserman. 

Remark 1.4. Consider the initial value system in n- dimensions (n > 3): 

L\u = |m| "-2 ; 
u{0) = Co; 
Vm(0) = ci. 

Theorem li.il applies to obtain some C^'" solution over a small neighborhood of 0, say, {x G : 
< R}. On the other hand, by a result of lOsf . if the solution exists in {x G M" : |a;| < R} 

2_ 2^ 

and Co > 0, then R < nu{0) "-2 = uCq ""^ . Consequently, i? — )■ as Cq +00. This does not 
contradict with Theorem \1.3\ apparently, since the solutions constructed in Theorem \1.3\ are of 
vanishing order 2m and hence Cq = 0. 

As a matter of fact, a large class of the systems fit into one or more of the three theorems. 
Especially, we have the following systems solvable. 

Example 1.5. For any p > 1 and any given R > 0, the system 

A'^u = ±\u\P 

has infinitely many C^"^'" non-radial solutions over {x G M" : |x| < R}, as a consequence of 
Theorem \1.3[ Here a = min{l — e,p — 1} with e any arbitrarily small positive number. Those 
solutions are necessarily smooth after a standard bootstrap argument. 

The following system has been well studied in the literature. 

Example 1.6. Let H E and H'{0) = 0. Consider the system 

Au = V{H{u)). 

According to Theorem \1.3\ for any R > 0, the above system has infinitely many non-radial 
solutions in C^'°'{{x G M" : |x| < R}) for any < a < 1. 

Indeed, a straight forward computation shows in the above example that a{u, Vm) = 
V{H{u)) = H'{u)Vu, Vpo{a{u,Vu)) = H"{u)Vu and Vpi(a(M, Vm)) = H'{u) and hence 
the system satisfies a G and a(0) = Va(0) = 0. By Theorem 11.31 for any i? > 0, there exist 
infinitely many solutions in the class of C^'"({x G M" : |x| < R\) and none of them is radially 
symmetric. 

One similarly has the solvability for the following m-th order Poisson type system. 
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Example 1.7. Let H e and H'{0) = 0. Then for any R>0, 

A™n = V{H{u, Vm, . . . , V^'"-^^)) 
has infinitely many non-radial smooth solutions in G M" : \x\ < R}) for any < a < 1. 

To see the solvability of the above example, a similar computation shows 

2m-2 

a{u, Vu,..., V''"-^m) = V{H{u, Vu,..., V'™-'m)) = ^ VjH{u, Vu,..., V^"'-\)V^+\ 

j=0 

where VjH is the derivative of H with respect to V-'m variable. Furthermore, 

2m-2 

Vp, {a{u, Vu,..., V''"-^m)) = J2 ^J^oH{u, Vu,..., V^'^-\)V^+\ 

j=0 

and for k > 1, 

2m-2 

Vp, {a{u, Vu,..., V^'^-^)) =Vp, ( Vj/f(u, Vu,..., V^'^-^u)V^+\) 

j=0 

= Yl ^j^kH{u, Vu,..., V'™-'u) V^+^n 

0<j,k<2m-2 

+ Vk-iH{u,Vu,...,V^"'-\). 

Hence a E and a(0) = Va(0) = 0. By Theorem 11.31 for any R > 0, there exist infinitely 
many solutions in the class of C^'"({x G M" : |x| < R}) and none of them is radially symmetric. 

The proof of the main theorems relies largely on a residue-type result for the fundamental 
solution r(-) of Laplacian in M", tt, > 3. This phenomenon is motivated by the Residue theorem 
for holomorphic functions in C In fact, the Cauchy integral formula specially indicates the 
integral of the Cauchy kernel with a degree k holomorphic polynomial over a simple closed 
curve in C is necessarily a holomorphic polynomial of the same degree. In this paper, we show 
the similar phenomenon also holds for r(-). Precisely speaking, denoting by Vk the space of 
polynomials of degree k restricted in {x G R" : |x| < R}, we have, making use of zonal spherical 
harmonics. 

Theorem 1.8. For any f ^Vk with k > 0, 




r{--y)f{y)dayeVk. 



Here day is the surface area element of {y G M" : \y\ = R}. 
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The above theorem plays an essential role in deriving the higher order derivatives for the 
Newtonian potential and the corresponding estimates via an induction process. Furthermore, 
as a side product, the residue formula allows to define the principle value of the higher order 
derivatives of the Newtonian potentials. 

The rest of the paper is outlined as follows. The notations for the function spaces and 
the corresponding norms are given in Section 2. In Section 3, we prove Theorem 11.81 As a 
consequence, the principle value of higher order derivatives of the Newtonian potential is well 
defined and computed in Section 4. As another application of the residue-type phenomenon 
for the fundamental solution, we derive the general higher order derivative formula and the 
corresponding estimates for the Newtonian potential in Section 5 and Section 6. Section 7 is 
devoted to the construction of the contraction map with the corresponding estimates necessary 
for the application of the fixed point theorem following the idea of |Pan2] . After a delicate 
chasing of the parameters, we show the main theorems hold in the last three sections. 

In Appendix A, a formula of the higher order derivative of the Newtonian potential over 
any general domain is derived making use of the same argument as in |GTj . Appendix B 
computes an interesting integral concerning the fundamental solution over the sphere making 
use of Gegenbauer polynomials. This provides a practical way to compute all the residue-type 
formulas for the fundamental solution. 

2 Notations 

Denote by B/j the open ball centered at the origin with radius i? in M", n > 3, and denote by 
OBr its boundary. Namely, B/j = {x G M" : |a;| < R} and dBji = {x E M" : |x| = R}. Here 
I ■ I is the standard Euclidean norm. We consider the following function spaces and norms over 
Bi^ following |Pan2j . 

Let C(Bji) be the set of continuous functions in B/j and C"(B^) the holder space in B/j 
with order a. For / e C"(B^), the norm of / is defined by 

II / |U: = || / II 

where 

II / lh= sup{|/(x) : X e Br}; 
HM] :=sup{^M£ip:x,x'GB^} 

if II / II a is finite. We note when || / ||q, is finite, the trivial extension of / onto B/j is then in 
C"(Bk). C"'(Bfi) is a Banach space under the norm || ■ ||q. 

For A; > 0, denote by C'^(B/j) the collection of all functions in B/j whose partial derivatives 
exist and is continuous up to order k. Denote by|| ■ \\ck the corresponding norm, where 

ll/llc. :=sup{|| D^f \\:\P\ = k} 
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if ll/llc* is finite. 

C^'°'{'Bji) is the subset of C^(B/j) whose fc-th order derivatives belong to C"(B/j). For 
any muhi-index (3 = . . . , /3„) with nonnegative entries, define \(3\ := YTj=iPj [3\ : = 
(3i\---(3n\- Given any / G C^(Br), we represent D^/ := 9f ^ ^2 ^ ■ ■ ■ " / with (9j- the partial 
derivative with respect to Xj variables. If / € C'^'"(Bij), we define the semi-norm 

||/||L^):=sup{|| D^f lU: = fc} 

if II / lli'^^ is finite. 

Of special interest, we introduce the subset of C^'"(B^) as follows. 

Co''"(B^) := {/ G ^^'"(B^) : D^/(0) = 0, |/3| < A; - 1}. 

The following lemmas play an important role in the rest of the paper. The proof can be 
found in |Pan2j and is omitted here. 

Lemma 2.1. [Pan2] If f E C'^'"(B^), then for any x,x' G Br and < a < 1, 

|/(x') -T,-(/)(x')| < C|x-xf+"( J2 Ha[D^f]). 

\fj,\=k 

Here T^{f){x') is the k-th order power series expansion of f at x. 
Lemma 2.2. [Pan2] If f E Co'''{Br), then for any I < k and < a < 1, 

||/||«<Ci?'=-' ||/||W. 
Remark 2.3. Lemma \2.S\ can be shown to hold for a = 0. 

As a consequence of Lemma [221 Co'"(Br) (0 < a < 1) becomes a Banach space under the 

II i\(k) 

norm || ■ ||a . 

Here and in the rest of the paper, we use C to represent any positive constant number 
dependent only on n, a and A^, where 0<a<l,n>3 and > 1. Especially, we point out 
that C is independent of R, which is later on a key parameter in the proof of the paper. 



3 Residue-type theorem for the fundamental solution in 

In complex analysis, the Residue theorem or the Cauchy integral formula states, for any holo- 
morphic function / in B^^ C C and any integer k > 1, z E B/j, 



f Ii^d^ = 2rrtf{z), 
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and hence 

J\^\=R '5^^ J\^\=Ry'i~^) 

where /*-*^^ is the fc-th derivative of / with respect to z. Recall the holomorphic kernel ^ is the 
Cauchy kernel for d operator in C and is also related to the first derivative of the fundamental 
solution in M^. As a special case, if / is a holomorphic polynomial of degree k in Bj^ C C, then 
for z e Br, 

[ D':^'{T^)fim = o- (3) 

When n > 3, there is no holomorphic kernel for d operator in general. 

On the other hand, the fundamental solution of Laplacian in M",n > 3 is 

T{x) = Cn 



\x 



n-2 ' 



Here c„ = „(2_^„)^ ? with Un the surface area of the unit sphere in M". With respect to the 
fundamental solution, the Newtonian potential of / in Br is defined by 

A/'(/)(x):= / T{x-y)f{y)dy 

for any integrable function / in B/j and for x G B/j. The Newtonian potential has attracted 
great attention in physics, and there have been many references concerning it for instance, |GT] 
and [NW] . 

The proof of Theorem 11.81 makes use of zonal spherical harmonics Zx^ and its reproducing 
property for spherical harmonics. See |SW] for the reference. In detail, let Hi be the set of all 
spherical harmonics of degree /, then for any f & Hi, 

fix) = [ Z^\y)f{y)day. 

Moreover, if / G H/^ with / ^ k, then 

0= / Zi'\y)f{y)day. 

On the other hand, denote by Vj^ the space of all homogeneous polynomials of degree k 
restricted in Br. For any / G V^, there exist -Pj's, some homogenous harmonic polynomials of 
degree j, such that 

f{x) = Pk{x) + \x\^Pk-2{x) + ■ ■ ■ + \x\^Pq{x), when k is even, (4) 

and 

f{x) = Pk{x) + \x\'^Pk-2{x) H h \x\^~'^Pi{x), when k is odd. (5) 
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Note PjldBi £ Hj. We now are in a position to prove the residue-type Theorem 11.81 for the 
fundamental solution F in M". 

Proof of Theorem 11.81 Without loss of generality, we assume / is a monomial of degree k. 
We also assume that R = 1. This is due to the following simple fact that for any / G Vj^ 



k ' 



/ r(x - y)f{y)day = R'+' [ r(| - y)f{y)day. 
Under the zonal spherical harmonics, we have, when x G Bi, 

r(x-y) = 5:c„,|-i^L_z|(^), 

where Cni = ^'+"'~^ . Letting y E (9Bi, the above expression for a; G Bi simplifies as 



T{x-y) = J2CnA4Z%{y). (6) 

1=0 

When k is odd, letting y G 9Bi and making use of ([5]), one has 

f{y) = Pk{y) + Pk-2{y) + ■ ■ ■ + Pi(y) (7) 

for some harmonic spherics Pj G Hj. Therefore, combining ([6]) and ([7]) together with the 
reproducing property of the zonal spherical harmonics, we have 

T{x - y)f{y)day = [ f V fp,(y) + Pk,,{y) + ■■■ + PM]day 



= a,fc|xrPfc(^) + Cn,k^2\xr'Pk-2{ri) + ■■■ + Cn,l\x\Pl{--) 

i-ly i/, lA,' 

= Cn,kPk{x) + Cn,k-2Pk-2{x) + " " " + C„,iPi(x) G Pfe. 

The case when k is even can be treated similarly and is omitted here. ■ 

We remark that, despite of the constructive proof of the Residue-type formula in Theorem 
ll.8l for the fundamental solution, the integral can actually be computed directly. See Appendix 
B for a computation of the formula when k = 1. The same method can practically be used for 
general k > 1. 

As an immediate consequence of Theorem II. 8[ we have the following corollary, analogous 
to ([3]) with respect to the Cauchy kernel in C. 

Corollary 3.1. For any f EVk and any multi-index /3 with > k + 1, 

I DlT{x-y)f{y)day = Q 

JdBn 

for X G B/j. 



8 



As another consequence of Theorem 11.81 we also have 



Theorem 3.2. For any f &Vk and any multi-index (3 with > k + 2, 

[ D^,T{x-y)f{y)dy = 0, (8) 

JBr\B,{z) 

when X G 'Qt{z) C B/j. Here 6^(2;) is the hall centered at z with radius e. 

Proof of Theorem 13.21 Write /3 = . . . , Without loss of generality, assume i? = 1, 
/3i > and / is a monomial of degree k. Moreover we write = — 1, . . . Hence 
applying Stokes' Theorem on V[x — y)f{y) over the domain B/j \ Be(z), one has 



! DlT{y-x)f{y)dy 

JBi\B,(^) 

-/ D^yV{y~x)d,f{y)dy+ f D^^Tiy - x) f{y)y,day (9) 

yi - zi 



D^;Tiy-x)fiy)^ \da^, 

9B,(2) \y - z' 



y 



Write / := /^^^ r(2/ - x)f{y)y,day and II := f^^^^^^ Dfr(y - x)f{yr^day. We show 
next that / = // in Bi and therefore 



/ D^yT{y - x)f{y)dy = - [ D^;V{y - x)d,f{y)dy. 

JBi\B,(2) JBi\B,(z) 

First note for J/, after a change of coordinates by letting y = z + er, 



(10) 



// =^-\^\ [ <r(^^ + T)fiz + er)nda^ 

JdBi e 

=e'-l^l / Dfr(^ + r)(/(er) + P,_i(r))rirfa. (11) 

JdBi e 

=e^-l^l+'= / D^jT{^ + r)/(r)nrfa.. 

Here Pfc_i(-) is some polynomial of degree k — 1. The last identity is due to the fact that / is 
a monomial together with an application of Corollary 13.11 onto Pa:-i(t)ti. 

When 1/3 1 > k + 3 and hence > k + 2, I and // are both zero due to the Theorem 11.81 
and we are done. When = k + 2, from ffTTj) we have 

// = / r(r + ^^)f{r)ndar. 



On the other hand, / by Theorem 11.81 is constant independent of x G Bi and therefore I = II 
when > k + 2 and hence (fTOj) holds. 

Now applying the induction process on (fTOl) . we get immediately, for a; G Bi, 

/ D^^T{y~x)f{y)dy = - [ D^^V{y - x)d,f{y)dy 

JBi\B,(2) JBi\B,(2) 

=CU) I D^yT{y-x)dy 

JBi\B,(z) 

=0. 

Here /i is some multi-index with |/i| > 2 and C{f) is some constant dependent only on / and 

4 Principle value of higher order derivatives of the New- 
tonian potential 

As a side product of Theorem 13. 2[ we show the principal value of the derivatives of the New- 
tonian potential exists. Denote by C^{W^) the set of smooth functions in R"- with compact 
supports. For any / G C^(M'^), recall the principle value of the derivatives of the Newtonian 
potential is defined as follows. 

Definition 4.1. p.v. J^^ L)^r(a; - y)f{y)dy := lim,_^o lRn\B,{x) D^'^ix - y)f{y)dy. 

Theorem 4.2. For any multi-index (3 with \f3\ = k > 2, f E , p.v. V^T^x — y)f{y)dy 
exists for all x G M^. Moreover, if suppf C B/j, then 

P-v. [ D^T{x-y)f{y)dy= [ D^T{x - y){f{y) - T^_,{f){y))dy, (12) 
jm" Jbr 

where T^{f){y) is the j-th order power series expansion of f at x. Moreover, the right hand 
integral of ( fl^) is independent of R. 

Proof of Theorem 14.21 We first show the independence of the right hand side of (JT2i) of R. 
Indeed, for any two numbers R' > R> with suppf C B/j C B/^/ and when x G suppf, 

T^.2U){y))dy+ I D^r{x-y){f{y)-T^_,{f){y))dy 
T^.2{f){y))dy - / D^T{x - y)T^_,{f){y)dy 
10 



/ D^Tix 


-y){f{y) 




[ D^T{x 

J Br 


-y){f{y) 


1 D^T{x 


-y){f{y) 



and hence by Theorem I3.2[ when x G suppf C B/j C B/j/, 

D^T{x-y){f{y)-T^_,{f){y))dy= [ D^T{x - y){f{y) ~ n_,{f){y))dy. 

\{ X ^ suppf , then 



/ D^V{x 


-y){f{y)- 


-n-2U){y))dy 


[ D^T{x 


-y)f{y)dy 




[ D^T{x - 

J Br 


-y)f{y)dy 




[ D^T{x ■ 
Jbr 


-y){f{y)- 


n_,U){y))dy. 



To prove ( fT2l) . we first note ( fT2|) is trivially true if x ^ suppf . 

When X G suppf , the right hand side of ( IT2l) is finite, due to the simple fact that when 
y e Bi(x), 

\D^V{x-y)\ < 

and 

1/(1/) -T,%(/)(i/)l<c(/)k- 1/1'-' 

with C{f) some constant dependent on /. 

On the other side, making use of Theorem I3.2[ one has 

Df'T{x-y)T^_,{f)iy)dy = 0. 



IBr\B,(x) 

As a combination of the above two facts, we obtain when x G suppf, 



p.v. [ D^T{x - y)f{y)dy = lim [ D^T{x - y)f{y)dy 

JM" JBr\B^{x) 

= lim / D^r{x-y){f{y)~T,^_,{f){y))dy 

^^"JBr\BAx) 



'Br\B,{x) 

D^T{x-y){f{y)-T^_,{f){y))dyl 
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5 Higher order derivatives of the Newtonian potential 



It is well known that the second order derivatives of the Newtonian potential in general do not 
exist due to the nonintegrability of the fundamental solution after differentiation more than 
once. However, when / is nice enough in the sense that / G C°(B/j), then N{f) G C^'"(B^). 
Especially, one has the following formula. 

Lemma 5.1. [Fr] Let f G C^CBr). Then for any x G ^r, 

d,d,M{f){x) = [ d.AT{x - y)U\y) - f{x))dy - ^/(a;). 
Moreover, for all f G C"(B/j), 

II W) lli'^<c||/|U, 

whenever || / ||o is finite. 

For higher order derivatives of the Newtonian potential, in the case when n = 2, |Panl] 
studied it from the point of view of complex analysis. There are few references in the literature 
for n > 3, though. In this section, we derive higher order derivatives of the Newtonian potential 
for functions in the appropriate spaces. 

Throughout the rest of the paper, unless otherwise indicated, we always regard derivatives 
inside the integration as derivatives with respect to y variables. For instance, inside an integral, 
d^Tix - y) := while d^Jix - y) := 

Definition 5.2. For a given multi-index (5 with = A; + 2, A; > 0, Hp : C^'"(B/j) — )• C(B/j) is 
defined as follows. 

MpU){x):= [ Dfr(x-y)(/(y)-T,-(/)(y))d|/, 
for f G C'^'"(Bij) and x G B/j, where Tj^{f){y) is the k-th order power series expansion of f at 

X. 

It is clear that the operator Affs is well defined over C'^'"(Bij). 

We next introduce the following notation for the convenience of the statement of the the- 
orem. Given any two multi-indices (3 = . . . and fi = (/ii, . . . ,yU„), we say /3 < /i if 
I3j < fij for each I < j < n and \(3\ < \fi\. Moreover, we define — (3 := {fii — ■ ■ ■ , fin — Pn) 
if /3 < fi. If in addition = + 1 with diD'^ = D^, we write fi — (3 = i . 

Definition 5.3. • Given a multi-index (3, we call j/?'--'^}^^]^ a continuously increasing 
nesting of length k for (3 if |/3^"'-'| = j for 1 < j < k and (3^^^ < (3^^'^^'^ < (3 for 
l<J<k-l. 
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• Given two multi-indices 7 and 7', we say 7' is the dual of 7 with respect to fi if D'^ = 

Making use of Theorem 11.81 together with Theorem IA.21 the foUowing theorem gives the 
formula for higher order derivatives of the Newtonian potential. 

Theorem 5.4. Let f E C'''°'{Br). Let (3 be a multi-index with \(3\ = k + 2 and a 
continuously increasing nesting of length k + 2 for (3. Then D^J\f{f){x) exists for x G B/j. 
Moreover, 



k+2 



j=2 |H=j-2 ^' 

Here 13^^^ is the dual of (3^^^ with respect to (3, and C{/3^^-^\ fi, ^^^^ - /3'-^-^^) is some constant 
dependent only on {(3^^~^\ n, (3^^^ — (3^^~^^). 

We point out, on the right hand side of ( |T3l) . the order |/i + of the derivative of / in 
the second term is always equal to k by definition. 

Proof of Theorem 15.41 = is given by Lemma [5.11 When A; > 0, for any multi-indices /3 
and /i with = + 1, one has by Corollary 13. H 

XB«(/3,/i,j)(^) := / D^^T{x-y)iy-xriy^day 

JdBn 



1 

R 
1 

R 



D^^r{x-y)y^y,day 



/ r{x~y)y'^y,day 

^l+H^^^ I T{^-y)y^y^day. 



According to Theorem 11.81 J^^^ r(^ — y)y^yjdijy is a polynomial of degree + 1 in x when 
\x\ < R and hence 

2:Bfl(/3,/i, = C{(3,n,j) 
with C{(3,fi,j) some constant dependent only on {/3,fj,,j). Therefore from Theorem IA.2I by 
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choosing Q = and Q' = B/j/ with R' > R > 0, one obtains 

D^Af{f){x) = I D^J{x-y){f{y)-n{f){y))dy 



k+2 



i=2 M=j-2 ^' 



Bp, 



D^^T{x-y){f{y)-T^{f){y))dy 



R' 

k+2 



J2 E ^^^^'''^'^'f^~^^'"'^^ z?^+^"^V(x) 



j=2 |/.|=i-2 

for any x G B/j and any R' > R. We then get N'{f) G C^+^(B/j). Moreover, for any /3 with 
1/3 1 < /c + 2 and for any x G B^^, after passing R! to i?, 

To simphfy the notation, we define an operator Tp by 
for any / G C^'"(Br). Then 7^ : C^'"(Br) ^ C:"(Br) and 

II r,(/) iu< c II / ||« . (14) 

Under this definition. Theorem 15.41 can be rewritten as 

D^uu) = mf)-r,U)- (15) 

for any / G C^'"(Br). 



6 Holder norm of D^M{f) 



In the derivation of the higher derivative formula of the Newtonian potential, the following 
operator shows up frequently and in itself is of an independent importance as well. 
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Definition 6.1. Given multi-indices (3 and /?' with \l3\ = k + 2, k > and = djD^' , a 
linear operator Sp : C''''"(Bj:j) — )• C(Bj:j) is defined as follows: 

Mf)ix) := / r(x - y)f{y)v,day 

for f G C'^'"(Bij) and x G ^r. Here day is the surface area element of OBr with the unit outer 
normal {ui, . . . , z/„). 

We point out Sf3{f) = T{-—y)f{y)vjday) is the counterpart in R*^ of the derivatives 

of the Cauchy integral for holomorphic functions in C For the convenience of computation, we 
shghtly modify the operator and define the following operator. 

Definition 6.2. Given multi-indices (5 and with \f3\ = k + 2, k > ^ and = djD^^' , a 
linear operator : C^^'^iY^R) — j- C(B/{) is defined as follows: 

Sp{f){x):= [ DPJr{x-y){f{y)~T^{f)iy)hday 

for f G C'''"'(Bji) and x G B/j. Here T^{f){y) is the k-th order Taylor series expansion of f at 
X, day is the surface area element of OBr with the unit outer normal (z/i, . . . , 

Note that due to Corollary O and = k + 1, Jg^ D^'T{x-y)T^{f){y)ujday ClS db function 

of X G Bij is a constant independent of R. As a result of this, i5/3(/) and S{f) are differed by 
a constant only dependent on Tj^{f), and especially, independent of R. 

We now are ready to prove the induction formula for the derivatives of the Newtonian 
potential. 

Lemma 6.3. Let f G C^'"(B^). let /3, (3' be two multi-index with \(3\ = k + 2 and = djD^' . 
We have 

D^UU) = DP'U{d,f) - SpU) - TpU)- 
Proof of Lemma 16. 3t Making use of Stokes' Theorem in (fT5|) and Corollary 13.11 we have for 
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Df'MU^x) =lim / Dfr(x - y){f{y) - T^{f){y))dy - r^(/)(x) 

= lim / D^Jr{x-y)d,{f{y)-T^{f){y))dy 

-lim / oJD^Jr{x-y){f{y)-T,^{f){y)))dy-rp{f){x) 
= lim / D^'r{x-y){d,f{y)-T^_,{d,f){y))dy 

- / D^Jr{x-y){f{y)-T^{f){y))u,day 

+ lim / r(x - - nUmWay - TpUm 

=D^W(9,/)(x) - - r,(/)(x). 

The last identity is because the third term is (9(e2-ra-fe-i+A;+a+n-i~) _ 0(e°). ■ 

The following lemma shows the operator Sj^ is a bounded operator from C'^'"(B/j) into 
C"(B^). 

Lemma 6.4. /ei (3 be a multi-index with \/3\ = k + 2. T/ie operator 5^ sends C'^'"(B/j) mto 
C"(Bh). Moreot;er, /or an?/ / G C'='"(Bij), 

IU< C7 II / ||('^), 

whenever \\ f \\^a^ is finite. 

In order to prove Lemma I6.4[ we need the following lemma. 
Lemma 6.5. For any x G Bi, < a < 1, 

/ 1 day<C(l-\x\r-\ 

Proof of Lemma 16.51 Assume x = (r, 0, . . . , 0) after rotation if necessary. One can assume 
in addition that r > ^. 

Choose spherical coordinates yi = cos 9i,y2 = sin 9i cos ^2, • • • , 1/n = sin 6*1 sin 6*2 ■■ ■ sin 6'„_i, 
where < < tt, < 6'j < 27r for 2 < i < n — 1 and denote by (9B"~^ the unit sphere in M"^^, 
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then we have, 



; day = / d6^ / daz 

Jo [(cos^i -r)2 + sin2^i]'^ Jdur' 

.c r- -""^ 

[1 — 2r cos + 



2 

[(1 -r)2 + 4rsm2|]V 
/ 7^ ^ d9+ / - — ^- rf^ 



=A + B. 
For A, since sin^ < 6 when > 0, 



For B, making use of sin^ > C9 when < < ^ and the assumption r > |, we get 

B<C de = C / sin--2 -d9<C (-r~'d9 < C(l-r)"-i+C < C(l-r)"-\ 

- yi_, (sinf)"- 2 - A-r^2^ - I ^ - K ) 

The lemma is thus concluded. I 

The following lemma in |NW] . Appendix 6.2a plays an essential role in the proof of Lemma 

Lemma 6.6. [NW] If z and z' are two points of the open unit disk in C, and 7 is the shorter 
segment of the circle through z and z' and orthogonal to the unit circle, then 



(1 — wwY ° 1 — a 
for < a < 1. 

Proof of Lemma [631 Write g{y) := f{y) - T^if){y). 
(i) The estimate for || ||. Indeed, by Lemma \2A\ 

I / Dl'r{x - y)g{y)u,day\ <C \\ f || « i?"" / \y 



--C II / ||(f) R-" [ \y- x\^-''+''day 

-cwfwi'^ I iy-fr-"^"rfa. 



<c II / Hi') 
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(ii) Given x^x' G Br, we estimate |iS/3(/)(x) — Sp{f){x')\. Assume without loss of generality 
that X, lie on the plane {y^ = ■ ■ ■ = Z/n = 0} and write x = Rz,x' = Rz' with z^z' G Bi. 
Then 

SpUm - S^UW) = [ (^f r(x -y)- D^'rix' - y))g{y)u,day 



--R-^ [ {D^jT{z -y)- D^Jviz' - y))g{Ry)u,day. 

JdBi 



Let 7(t) = (7i()f:), 72(t), 0, . . . , 0) : [0, 1] — )■ {y^ = ■ ■ ■ = = 0} = C be a parametrization of 
the shorter segment of the circle through z and z' and orthogonal to the unit circle in C with 
7(0) = ^,7(1) = z. We then have 



Sp{f){x) - Sp{f){x') =R'' [ [ l(D^'r(7(t) - y))dtg{Ry)v,day 

JdBi Jo O.'^ 

=R-^ ['f^^'^m [ r(7(t) - y))g{Ry)u,day. 

Jo JdBi 

Making use of Corollary 13. we have for any < t < 1, 

/ (9,,D^'r(7(t) - y))g{Ry)u,day = [ {d,,D^;r{^{t) - y)) {g{Ry) - T^^^'\g){Ry))u,day, 

JdBi JdBi 

where T^"'^^\g){y) is the k-th order power series expansion of g at R'y{t). Furthermore, by 



Lemma I2.H 

\g{Ry) - T^'^^'\g){Ry)\ < C\Ry - i?7(t)l'+" E ^-[^"^] 

/i|=fc 

= -7(t)i'+" Ha[D^n- 

iMi=fc 

Therefore, 

bit) - - 7(t)|'+"rfcT, 



|M|=fc fc=l -^^Bi 

\u\=k k=l -^dBi 



1 2 

^{t)-yr^'^day 
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Applying Lemma [631 to /^^^ |7(t) — y\ "'^°'day in the last expression, 



Hence by Lemma 16.61 



-\iitwr 

\dw\ 



\Sp{f){x) - Sp{f){x')\ <C{J2 H4D^f])R'^\z -z'r = C{J2 H^[D^f]) \x - x'^- 

|/i|=fc \tA—^ 

Namely, 

\t,\=k 

We finally have shown, combining (i) and (ii), 

II SpU) IU< c II / ||(f) .1 

Remark 6.7. In the proof of Lemma \6.4\ ii) when carrying out the holder norm for Sp, a 
natural choice of 'y would usually be the segment connecting z and z. However, the estimate in 
Lemma \6. 6\ actually fails if is chosen to be the segment instead of the geodesic as in INWf . 

Applying Lemma 16.31 and Lemma 16.41 inductively, one eventually obtains the following for- 
mula. 

Theorem 6.8. Given a multi-index (5 with \(3\ = k + 2, let {(i^'^^} be a continuously increasing 
nesting for (3 of length k + 2 and (3^^^ be the dual of 13^^^ with respect to (3 for 2 < j < k + 2. 
Then for any f e ^='"(6^), 

D^Mif) = D^'^WiD^'^^' f) - X^5,o)(D^"''/) - r,(/), 

i=3 

in Bij. Moreover, for any f G C'''"{'Bji), 

II m) iiL'^'^< c II / 1|(^'), 

and consequently, for any m G 

II Ar™(/) ||i'=+2^)< C II / ||(f), 

whenever \\ f \\a^ is finite. 
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Proof of Lemma 16.81 By Lemma 16. 3[ 



k-\-2 

i=3 

Hence from the above identity, for any / G C'^'"(Br) as long as || / is finite, 

\\Afif)\\'i+^y.= sup \\D^Af{f)\\^ 

|/3|=fc+2 

k+2 

< C sup [ II D^''W{D^'''' f) lU II S^u,{D^'''' f) lU + II Mf) lU ]• 

Since \j3^^^\ = j and = k + 2 — j from definition, by Lemma [5 .It Lemma [6.41 and ( fT4l) . we 

get 

fcH-2 

II W) 11^'^ < C sup [ II MiD^^^'f) IlL^) + X: II D^^'f ^t'^ + II / lli'^ ] 

< C sup [ II D^'^'f lU + II / ||(f+^-^-+^-^) + II / lli^) ] 

|/3|=fc+2 

< II / IlL'^ • 

Finally, applying induction in the above expression, 

II -^"(Z) 1^+2"^) =11 J\f(J\f^~^(f)) II ^=+2™) 

< C II J^f^^^^f) ||(^+2m-2) 

< • • • 

< II / lli'^ •■ 

7 Construction of the contraction map 

In this section, we construct a contraction map from the system ([T]). Assume a G C^. For any 
vector-valued function / G (Co'"'"(Bij))^, introduce u^^^f) := {u^\f), . . . ,u^j^\f)) with 

^f\f){^)= I T{x-y)a,{yJ{y),Vf{y),...,V'"'f{y))dy 
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for 1 < j < A^. According to Theorem EH u'^^^f) E (^^'"(Bij))^ and 



IIL'^<C|| a,(-, /,..., V^'"/) lU 



We define inductively = {uf \f ),..., J^\f)) for 1 < / < m as follows. For each 

1 < j < and a; G Br, 

cof\f){x) ■.= ^^{cof~'\mx). 

Note that, in terms of the Newtonian potential, 

uf\f)=^^'{a,{.,f,...,v''-f)). 

Therefore, by Theorem EH co^^\f) E (^^''"(B^))^ and 

ll^f(/)lli''^<l|a.(-,/,---,V^™/)|U. (16) 

We also define 9{f) := (6'i(/), . . . ,6'iv(/)) from uj^"^\f) by truncating degree less than 2m 
terms and part of the degree 2m terms in its power series expansion at 0. Precisely speaking, 
for 1 < j < A^ and x E Br, 

om^) = - T,_,^-)(/))(x) - y: ^Mll^M^/^, (17) 

/3eA ^' 

where T2m-i{^j^\f)) is the (2m — l)-th power series expansion of u}j^\f) at 0, A = {/3 : \/3\ = 
2m, and at least one of (3j is odd for 1 < j < n}. 

From the construction, it is immediate to see that for any / E (Co™''"(Br))^, A'^9{f){x) = 
a(a;,/(x),V/(x),...,V2™/(a;)) when x E Br. Moreover, u^'^^f) E (C2'"'"(Bh))^ and so 
elf) G JCo"''"(B^))^. We note that, because of ([16]) and j^, e{f) is automatically in 
(Cq™''"(B/j))^ after a trivial extension onto Br if f E {Cq"^'°'{Br))^ . 

Recall {Cq^'°'{Br), || ■ \\a^) is a Banach space. We now have constructed an operator between 
two Banach spaces as follows. 

e : (Co^'"'"(B^))^ -> (Co'™'°(Bk))^ 

with the corresponding norm 

II f ||(2m)_ II f ||(2m) 

" ^ "° ~ l<i<^ " "° ■ 

The ball of radius 7 in Cq"^'°'{Br)^ is denoted by 

B{R,^) := {/ E Co^'"'"(B^))^ :|| / ||(2™)< 7}. 
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On the other hand, recall a function u G C^'^ is called fc-harmonic if A'^u = 0. Given 
h = {hi, . . . ,h]\r) with hj any homogeneous m-harmonic polynomial of degree 2m and for any 
/ G (Co''"'°(B^))^, consider 

9,{f) = h + 9{f). 

Then e,{f) G (^'"(B^))^, A-^,(/)(a;) = A™e(/)(a;) = a(x, /(x), V/(x), . . . , V2-/(x)) in 
B/j while the 2m jets D'^9h{f){0) with /3 G A coincide with those of the given /i. 

We will seek the solutions to ([1]) by making use of the fixed point theorem. Indeed, we first 
show there exists 7 > and R > 0, such that 9 : B{R, 7) — )■ B{R, |) and 6^ is a contraction 
map. We then pick some nontrivial h as above with h G B{R, ^) and consider the corresponding 
operator 9h- Consequently, 9h '■ B^R,^) — )■ B{R,'y) and is a contraction map. As an application 
of the fixed point theorem, there exists some u G (Cq'"'"(Bj:j))^ such that 9h{u) = u. This u 
apparently satisfies A'^m = A"^9h{u) = a{-,u, Vm, . . . , V^^u) in Bjj and hence is a solution to 
([T]) over Br. 

Remark 7.1. The construction of 9^ guarantees the solution u obtained from the fixed point 
theorem is not a trivial solution. More precisely, u is of vanishing order 2m. 

We divide our proof into two steps. In each of the steps, we need to utilize Theorem 16.81 
and then the estimates in |Pan2] . 

7.1 Estimate of || e{f) - e{g) H?""^ 

First, we note from ( fT71) that for 1 < J < A^, for any f,g& B{R, 7), 

II dAf) - 0A9) lir^< II -^{f) - u^fHg) 11^^ + II v^-(c.f )(/) - II 

<\\jr\f)-^''r\9) iir^ 

= 11 M{uf'-'\f))-M{Jr'\9)) iii'™^ 

= II Ar'"(a,(-, /(■), V/(-), . . . , V'^^fi-)) - a,{;g{-), Vgi-), . . . , V'^^gi-))) HL^'") 
Making use of Theorem 16.81 into the above expression, we have then 

II OAf) - d^g) lli'™^<C^ II /(■), V/(-), . . . , V2"/(-)) - a,(-, (?(■), Vg{-), . . . , V^'-gi-)) |U ■ 

(18) 

We next proceed with the same derivation of the estimate as in section 4.1 of |Pan2] without 
proof. Note that due to Lemma E^l when / G -B(i?,7), then || II < CR^'^-^-f for < j < 
2m. Therefore as a vector-valued function of {p-i,Po,Pi, ■ ■ ■ ,P2m), cl takes value in E := {p-i G 
^R)Pj ^ Bc^2m-:i^, < j < 2m} when u G B{R, 7). 
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Denote by Aj := sup^ I Vp^.a|, Qj := sup { ^^^^^^^z^^^^^ : x,x' e E} and Lj := sup£;(V^^.p2^a) 
with — 1 < j < 2m. Therefore, for — 1 < j < 2m, 

< C II Vpj.a ||c(£;)< C* || a \\c^,<^(e), 

Qj < C* II a ||ci."{£;), (19) 
Lj < C II Vp2^a \\c^{E)< C II a ||c2(£;) . 

Here || a ||ci,a(£;) = || a ||ci + sup { ^^^^^|^f^^^ : G 

Lemma 7.2. [Pan2] For any f,g E B{R,'j), if a e C^, 

II a{; /(■), V/(-), . . . , V'^^fi-)) - a{; g{-), Vg{-), V'^gi-)) |U< S{R, i) \\ f - g |U, 

where 

2m 

6{R, l) = CY, R^"""' {^3 + ^"(1 + R'l" + l)Qj + iLj). (20) 
i=o 

Moreover, if a is independent of p2m, then when a G C"^'", 

2m- 1 

5(i?, 7) = C ^ (A,- + i?"(l + /2"7" + 7)Q,) . (21) 

i=o 

We then have obtained from ( fTSl) . by using Lemma [7.21 that 

II W-%)llf"^<Wr)||/-^?IU, (22) 
with 5(i?,7) given in (EDI) or (EI]). 

7.2 Estimate of || e{f) HJ?""^ 
Similarly, for / G i3(/?,7), 1 < j < A^, 

II %(/) lli''"^< II ^r(/) lli'"^^ +|V^-(a;f )(/))| 

< II ^t\f) lli'"^ 

= l|A/''"(%(-,/(-),v/(-),...,v^-/(-))) llL^"^) 

<C II «,•(-,/(■), V/(-),...,V2'"/(-)) lU- 
According to the estimate in section 4.2 of jPan2j . 



(23) 
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Lemma 7.3. [Pan2] For any f E B{R,'j), if a e C^, 

II a(-,/(-),V/(-),...,V2'"/(-)) |U< r/(i?,r), 

where 

r]{R, 7) = |a(0) I + C{R{A^^ + (l + i?"7" + l)Q-i + 7^-i) + l)) (24) 

with 7) g'zven in Ij^EE)- 



Moreover, if a is independent of p2m, then when a G 0^'°", 

r/(i?,7) = |a(0)| + C{R{A^i + + + l)Q~i) + 7<5(^,7))- (25) 
wt/i 5{R, 7) g'Zf en m (fH])- 

Combining Lemma 17.31 and fl23|) , we have 

II ^(/) ||L''"^<r^(i?,r), (26) 

with ri{R,j) given in or in fl25|) . 



8 Proof of Theorem 11.2 



We now prove a shghtly more general result than the main theorems following jPan2j . 

Theorem 8.1. Let a e and a(0) = 0. There is a constant 6{< 1) depending only on n, N 
and a, such that when 



|Vp,„a(0)| + |V^,„,,„a(0)|<5, 

the system ([!]) has infinitely many solutions in C^*"'" of vanishing order 2m at the origin in 
some small neighborhood. 

Proof of Theorem 18.11 Our goal is to show 6 sends B{R,'y) into B{R, ^) for some positive 
R and 7 and is a contraction map between B{R, 7). In other words, we show there exist 7 > 
and R> such that for any f,gE B{R, 7), 

II 0{f) - e{g) ||i2™)< c II / - ^7 IlL'") with c < 1 

and 

II Oif) llf"^< I- 

From ( I22I) and ( l26l) . it boils down to show there exist 7 > and R> such that 

S{R,-f) < c < 1 
(J, ^27) 
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Denote by r := |Vp2„fl(0)| + iVp^^p^^'^lO)!' ^li^) represent a constant converging to 
as -R — )■ for each fixed 7, and use e{R + 7) to represent a constant converging to as both 
R and 7 go to 0. Then by continuity of a, 

< r + e(i? + 7), Q2m<CT + e{R + ^), < r + e(i? + 7). 

fl20|) and fl24|) can hence be written as 

7) = C„r(l + 7) + e^(i?) + e(i? + 7), (28) 
r]iR,^)=Cal6iR,^) + e,iR). (29) 

with Ca dependent on || a ||c2(e)- 

First, for each 7, choose Rq such that e^{R) < j when R < Rq in Then we will choose 
7 and i? small enough so (5(i?, 7) < c := minj^^r^, i} < 1. Indeed, by choosing 7(< 1) and R 
small, we can make e^{R) + e{R + 7) < | in (128|) and hence 

S{R,j)<2CaT+^. 

When r < we thus have (!27jl holds. 

Now recall A = {/3 : = 2m, and at least one of (3j is odd for 1 < j < n}. For R 
and 7 chosen as above. Pick h{x) = bx^ with /3 G A, and make 6 > small enough such 
that II h ||q'"''< ^ and hence h G 13{R, |). Consider the operator 9h{f) := h + 6{f)- Then 
6h '■ S{R, 7) — )■ i3(-R, 7) forms a contraction map from the construction. By fixed point theorem 
for Banach spaces, there is some u G B{R, 7) such that 9h{u) = u. u thus solves the system ([1]) 
in the class C^*"'"" and is of vanishing order 2m by the construction. ■ 

Remark 8.2. None of the solutions constructed in the proof of Theorem \8.1\ is radially sym- 
metric, i.e., none of them is obtained by reducing the system ^ possibly into an ODE system 
with respect to the radial variable r = \x\ only. Indeed, if the solution u{x) = u{r) G Cg™'", then 
near 0, u{r) = er^™' + odrp"^) for some constant e. In particular, D'^u^O) = for all /3 E A. 
This apparently can not happen because from the construction, h = bx^° with some Pq E A and 
D*m(0) = D'3«/i(0) ^ 0. 

Proof of Theorem 11.21 Theorem 11.21 is a consequence of Theorem 18.11 and Remark 18. 2[ I 



9 Proof of Theorem 1.1 



When Cj = 0, < j < 2m — 1 and a is independent of p2m, ^2m, Q2m and Lj{—1 < j < 2m) 



are all and so (l22l) and (1261) becomes 



5(i?,7)<e,(i?), 
r/(i?,7)<|a(0)| + e^(i?). 
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Here we only need C^'" regularity for a from the estimates (12T]) and (1251) . Now we choose some 
positive 7o so that 70 > 4|a(0)|. Consequently, we choose R sufficiently small so e-yg{R) < c := 
min{|,^} < 1. Hence 

S{R,jo) <c<l; 



Applying the same strategy as in the proof of Theorem 18 -H we can find a solution u G B{R, 70) 
to the ODE system ([2]) which is not radially symmetric. 

For general given C/3's with multi-indices /3, we write T2m-i{x) '■= Yl^=o^ ■ Consider 
the new system 

N^u{x) = a{x, u + T2„_i(x), V(m + T2^_i(x)), . . . , V^'^-\u + T2m-i{x))- 
L)/5m(0) = 0, 0<|/3|<2m-l. 

This is a system with all the initial values equal to 0. We then obtain some solution u in the 
class of C^™"'" in some small neighborhood of 0. Then u = u + T2m-i solves the system ([2]) in 
the class of C^™'"" in some small neighborhood of 0. Apparently, the solution obtained in this 
way is of vanishing order at most 2m. Moreover, u is not radially symmetric since u is not. I 



10 Proof of Theorem 11.3 



Since a is independent of x and a(0) = 0, v4_i, Q-i and L_i are and hence in 

r]{R.l)<Cal5{R,^). 



In order to prove Theorem 11.31 we need to show for any fixed i? > 0, there exists some 70 > 
such that 

W7o)<i; 



r/(i?, 7o) < y. 



which is equivalent to showing 



W7o)<c:=min{i,^}<l. (30) 

Indeed, since Va(0) = 0, we have a G Cq'^{E) and hence by Lemma for < j < 2m, 
A, < C II Vp^a ||ci(£;) R^^^-'i <C\\a \\c2(e) R^^^-'l- 
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On the other hand, we also have by definition, for < j < 2m, 




where 6/2(7) represents some function converging to as 7 goes to for each fixed R > 0. fl30l) 
is thus true and the proof of Theorem 11.31 is complete. I 

Appendix 

A Higher order derivatives of the Newtonian potential 

We derive the higher order derivatives of the Newtonian potential following |GT] . Let Q C M" 
be bounded. 

Definition A.l. Given two multi-indices /3 and fi, and j with 1 < j < n, we define for x G Q, 



where day is the surface area element of dQ with the unit outer normal (z/i, . . . , z/„). 

It is clear to see that XQ(/3,/i,j) G C°°(i7). 

Theorem A. 2. Let f3 be a multi-index with \f3\ = k + 2. Let {P^^^} he a continuously increasing 
nesting of length k + 2 for (3 and let (3^^^ be the dual of f3^^'> with respect to f3 for I < j < k + 2. 
Then given a bounded and locally C^'" function f in Q and for any x G fi, 



Here Q' D Q are such that Stokes' Theorem holds on Q' and f extends to vanish on Q' \ Q. 
Proof of Theorem IA.2[ The theorem is proved by induction on k. When k = 0, the theorem 



reduces to the case in |GTj . Assume ( 1311) is true for k = ko > 0, i.e., for any / G C'^°'", and any 
(3 with \(3\ = ko + 2, 





(31) 




(32) 



j=2 M=j-2 
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We want to show it is true ior k — ko + 1. Namely, for any /3 with \/3\ — ko + 3 and / e C'^o+i,"^ 



D^j\f{f){^) = / ^fr(x - y){f{y) - Ti^Af)i.y))dy 

Jo.' 



fco+3 DLL f ( \ ("^"^) 

3=2 \^l■\=3-2 ^' 



Without loss of generality, assume = diD'^^'"^^^ with = A;o + 2. Let 

ve{x) = I Z)f °"^^r(x - y)r^,{x - y) {f{y) - Tl{f){y))dy 
Jii' 



J=2 |M|=j-2 ^ 

where 7]^(x—y) = ry( ^^~^^ ) with some smooth increasing function such that r]{t) = when t < 1 

and 77(i) = 1 when i > 2. Here we choose e < When e ^ 0, v^{x) D^''''°^^''j\f{f){x) 

for all X e Q by induction. 
Now consider 

dMx) = - I d,{Df'^^''T{x-y)ve{x-y)){f{y)-T^^{f){y))dy 
+ f Df''^''T{x-yMx-y)d,,{f{y)-T^^{f){y))dy 

i=2 M=j-2 ^' 

Here A :^ - J^, d^iDf^'^'^Fix - y^x - y)) (/(y) - T^,{f){y))dy and B := J^, i^f''=°"^^r(x - 
y)77e(x - y)d^i{f{y) - T§^{f){y))dy. We will show as e 0, for all x en, 



A + B^ [ D^,r{x-y){f{y)-T^^^,{f){y))dy 
Jn' 



|/i|=feo + l ^' 
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3^ thus gives for x G f2, 

dM^) ^ / ^fr(x - y){f{y) - Ti^,um)dy 



i=2 |A«|=i-2 ^' 



and hence ([33]) is concluded. 
For A, 



A = - [ 5i(l^r"'r(x-y)r/,(a;-y))(/(y)-T,^,^i(/)(y))rf2/ 
Applying Stokes' Theorem to the second term of the above expression, we then have 



A = - d,{Df"^''v{x-y)rt,{x-y)){f{y)-Tl^,Um)dy 
- E RUM [ Dr^'^r{x-yMx-y){y-xru,day 

+ E RHIj^ f Df^'''Tix-yMx-y)My-xrdy. 



|M|=fco+l 

On the other hand, 

B = - [ Df'^''v{x-y)r^,{x-y)d,,{Tl{f){y))dy. 
Jci' 

Therefore 

A + B = - [ 9i(Z^f°"^'r(a:-y)r^,(x-y))(/(y)-T,^,,^,(/)(y))rfy 
Jn' 

- E I Df'^''T{x-y)^Xx-y){y-xYv,day 

+ [ Df'^''T{x-y)r^,{x~y)[ ^^d,{y - xT - d^,{T^^{f){y))]dy 

Jn' I I , , 1 Z^' 



--I + II + III. 
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As e 0, for x e ^l, 



[ D^,T{x-y){f{y)-Tl^,{f){y))dy 
Jn' 

E f Dt'''\{x-y){y-xY.,da, 

, . , , Jan' 



J y,yy ^.^.y ^gg^ 

/i|=fco+l 

D^'f{x 



|At|=fco+l 

For ///, notice T^; (/)(») = E|,|<i, ""-"t!"-"' . so 

\^^\<kQ ImI<A;o 
On the other hand, one observes 

diD>^f{x){y-xY ^ D^fix)^, , d,D^^fix){y-xr 



|/x|<fco ^ |/x|=fco + l ^ 

Hence 



^ ^^9i(y-xr-9.,(T,^„(/)(y))=0, 

|/i|=fco+l ^' 

and 

/// = 0. (37) 
Combining ([Ml) and dST]), ([35D thus holds. ■ 

B Computation of ^^^(O, 0, 1) 

We will compute Xb^(0, 0, l){x) := /^^^ V{x — y)viday for x G Bi. 

Write X = f/ ■ [a, 0, . . . , 0]*, where f/ = (Mij)i<i,j<n is some unitary matrix and a = \x\, and 
then make a change of coordinates by letting y = U ■ y m the expression of Xbj(0,0, 1). We 
then get 

Xb, (0,0,1) := / , -^, -2 
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Write y back by y allowing an abuse of notation, then 

Jij,(0,0, 1) :=uii / „ 2 *^0"^ 

JdB, ^(a - yif + yl + --- + yl 

Since .. —^ is odd with respect to yj when j > 2, 

A/(a-yi)2+j/2+-+j,2 

rrfcT,, = 



IdB, ^{a - y^f + yl + ... + y2^- 
when j > 2 and hence 

Xb,(0,0, 1) := Mil / „ 2 ^^v 

JdB, ^(a-y^y + yl + '-' + yl 

Next, rewrite the above integral in terms of the spherical coordinates, then we obtain 

_ , , sintcos"~^t , 

Xbi(0,0, 1) =w„„iMii / „ — 

J-^ v/(a-sint)2 + cos2t 

=w„_iuii / —^du 

J-i (1 -2aM + a2) — 

Here Un-i is the surface area of the unit sphere in R"~^. 

In order to compute fl38l) . we need to use Gegenbauer polynomials. Recall for each fixed p, 
the Gegenbauer polynomials are {Cn\x)} in [—1,1] C M satisfying 



^ oo 
^ ' n=0 



in (—1, 1). In particular, 

Ci,'\x) = 1, 
C[''\x) = 2px, 

Ci^\x) = ^[2x{n + p - l)dA{x) - + 2p - 2)CiU^)]. 

Moreover, {c!{'\x)} are orthogonal polynomials on the interval [-1,1] with respect to the weight 
function (1 — x'^Y~'2 . In other words, 

1 

CiP\x)Clll\x){l-xy-'^dx = {)^m^n 

-1 
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Letting p = then — ^j^rr = X]^o^"''''('")'^" ~ • can hence be 

(1— 2au+a2) 2 " 

written as 

■^"1 n=0 ^ 

J -I 

a 7r2i-2pr(l + 2p) 

-""-"-2;^ (1 + P)r(p)^ 

47r2 



-Xi. 



Remark B.l. Making use of the same approach as the above in addition to observing some 
symmetry of the integrand over the sphere, one can practically compute XB^(/3,/i,j) for all 
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